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The recent observation of the quantum state of the neutron bouncing freely under gravity aUows 
some novel experiments. A method of purifying the ground state is given, and possible applications 
to the measurement of the electric dipole moment of the neutron and the short distance behaviour 
of gravity are discussed. 
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I. INTRODUCTION: 

In an elegant experiment, Nesvizhevsky et al Q 
have observed the quantum states of a neutron bounc- 
ing freely under gravity. Their experiment shows that 
the state exists with approximately the predicted prop- 
erties by observing its extinction by lowering a neutron 
absorber. In this paper, we discuss a method for puri- 
fying the ground state and then discuss experiments to 
exploit it. 

Ultracold neutron (UCN) technology Q has been 
around for many years, and it is exploited in the lat- 
est EDM experiments This relies on the observation 
that (for example) beryllium metal totally reflects neu- 
trons at a temperature of below 100 K. This implies the 
existence of the quantum state; that of a neutron bounc- 
ing freely under gravity. Although the idea is a fairly 
obvious one, and the theory is very old |^, Hf, it appears 
that no attempts were made to observe it until recently. 
It seems possible that if the ground state of the system 
could be isolated, there are several experiments which 
could be performed 

In the next section we show that a simple device based 
on a classical centrifuge could be used to purify the 
ground state. We then show that the lifetime of the state 
is expected to be long. We then consider two possible 
applications to fundamental physics. Firstly, it is well- 
known that the breakdown of CP invariance and hence 
T invariance in fundamental interactions implies the ex- 
istence of an electric dipole moment (EDM) for the neu- 
tron. Predictions of the magnitude vary from 10~^^ e-cm 
down to 10^"^^ e-cm, while the current limit is 0.5 x 10~^^ 
e cm. Obviously it is of considerable interest to lower 
this, since it restricts possible extensions of the standard 
model considerably. Finally we investigate the possibility 
of using the system to test for extra gravitational inter- 
actions at very short range. 



II. THE BOUNCING NEUTRON 

In this section we review the theory for the one- 
dimensional case. The bounce eigenfunctions Zn{z) are 
defined by the Schrodinger equation which describes the 
neutron under gravity is 



2m dz^ 



Z (z) + (ctz - S) Z (z) = , cr = mg (1) 



This can be converted to a dimensionless form via the 
substitution 



,'2mcr 

y = pz-y„, (3^ [ 
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(2) 



The boundary condition for a totally reflecting " ground" 
is that the wave-function vanishes at z = 0. This gives 



with 



— Z„ {y) + {y- yn)Zn {y) = 



Zn (-yn) = 



(3) 



(4) 



As is well known, this is the equation for the Airy func- 
tion, and hence t/„ is the n-th zero of the function. This 
gives rise to an energy En and a "bounce height" z„ 



r, yn J Zn 

2m fi 
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For the neutron, the ground state parameters are Eq = 
lAOlpeV and zq = 13.7 ii (these results differ slightly 
from those quoted in 0]. The latter number implies 
that one could build a macroscopic apparatus sensitive 
to the parameters of a quantum state. The usual method 
of looking for bound states, by measuring the spectrum 
of emitted photons, is not available. The frequency for 
the 1-photon transition between the first excited and the 
ground state is 254 Hz 12J and one can estimate the ra- 
diative lifetime of the first excited state to be much larger 
than the Hubble time! 
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FIG. 1: The neutron centrifuge. A shows a classicaUy con- 
fined orbit, B shows a higher energy one that escapes. 

III. THE NEUTRON CENTRIFUGE 

The basic concept is shown in Fig ^UCN's are injected 
into the centre of the apparatus. The height of the step h 
is adjusted so that the lowest energy bounce states will be 
trapped in the middle of the apparatus while the higher 
energy ones will bounce over the step and will be cen- 
trifuged out. It is obvious that this will work classically: 
it is less obvious that it will work quantum-mechanically. 
In order to investigate this, we need to solve the 3-D 
Schrodinger equation which describes the neutron: 

V^* (f) + ((7z-E)'^{r) = , a = mg (6) 

2 m 

In the case where the step-height h = 0, the wave func- 
tion 5'm(r, z, (j>) can be separated to give 

^'fein(r,z,0) -/zfe(r)Z„(z)e^"^ (7) 

We now want to use these as the basis functions for 
the problem with a step. The step is defined in general 



via a function H{r): in this case 

H (r) = r < i?i H{r)=hr>Ri (8) 
We can then define the new wave function as 

$™(r,2;,0) = ^a;"i„(r,z)*fein(r,2;,0) (9) 

The solution of this problem is shown in the appendix. 

To simulate the passage of a large number of neutrons 
through the centrifuge, we first choose a random super- 
position of the <&„i(?', z). 

eo(r,z)=^6^„a>„(r,z) (10) 

After passing once through the centrifuge, all compo- 
nents of the wave function lying beyond r = Rg are elim- 
inated, so we find a new wave function 

e,(r,2) = ^6:,,$„(r,z) (11) 

where 

6:„ = ^d™™'6L"' (12) 

and 

drum' = J ^rn{r, z)^m'ir, z)9 {Rs - r) rdrdz (13) 

Note that the dram' only depend on the eigenfunctions, 
so it acts as a transfer matrix. The phases of the V^n ^"^^ 
randomized for each pass round the centrifuge. Strictly 
speaking, this is not adequate: instead we should form 
wave-packets and project out a superposition of the ra- 
dial and bounce states. However, the process we have 
defined here is essentially equivalent. 

To obtain actual results, we must choose a set of pa- 
rameters. We use Ri = lO/i, R2 = 20^, i?3 = 35/i, Rsut = 
30/i and h = VJ fx: these are fairly arbitrary, with the ex- 
ception of h which must lie between the bounce heights 
for the ground and first excited states. The set of basis 
functions includes 12 bounce states and 6 radial states, 
for a total of 72: again this is not a critical parameter. 

The process is very efficient: in fig |5] we show the 
initial probability Pa{r,z) = |6o(7', z)]^ In in fig 13 we 
show the corresponding probability Pioo(7', z) after 100 
passes. Obviously the wave function is very strongly con- 
centrated in the central region, and corresponds on the 
lowest bounce states (n = 1). It is not straightforward to 
quantify this further, since our actual wavefunctions are 
a mixture of various bounce and radial states. However, 
it is possible to find the mean height after each passage 
through the apparatus, and compare this to the mean 
height of the lowest bounce state, which turns out to be 
9.6 /i. This is shown in fig ^ where we plot the mean 
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FIG. 2: initial probability distribution 




FIG. 3: Probability distribution after 100 passes 




FIG. 4: Mean Height of state as a function of the number of 
bounces 

height as a function of the number of passes. These re- 
sults are for 1=1: they are in fact more dramatic for 
larger 1, as intuition would suggest. 

We conclude that this entirely passive device will per- 
form a very effective job of separating the lowest bounce 
state. Two points should be noted, however. In practice 
we will not know if there is a neutron in the device. Sec- 
ondly the step effectively mixes bounce and radial states 
and their energies, so that (e.g) each bounce state would 
be spread over a number of radial states. Careful adjust- 
ment of the radial parameters can minimize this. 

IV. FINITE PENETRATION EFFECTS 

We have assumed above that the neutron is perfectly 
reflected. In fact this is not correct: the potential barrier 
at the metal has a finite height, given by the Fermi quasi- 
potential (3), which allows the wave- function to pene- 
trate into the surfacep^. This has two effects: firstly the 
calculated energies and bounce heights will not be cor- 
rect. Secondly the penetration of the neutron allows it to 
be either absorbed or scattered inelastically by phonons: 
obviously either of these would destroy the state. Hence 
even the ground state would be expected to have a finite 
lifetime: if the lifetime is shorter than the neutron life- 
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time, it would seriously affect any potential experiments, 
so it is important to estimate it. 

Provided that the penetration depth is small, we can 
estimate these effects by perturbation theory. The po- 
tential inside the metal is 



SV = 6Vr + iSVi 



(14) 



where 



SVr = 252neV, SVi = 1.26peV (15) 

are the values appropriate for Be other choices of 
material do not significantly affect the argument. The 
wave function for the n'th excited state inside the surface 
is given by 



Zn (z) = A„e""" (z < 0) 



with 



En) 



(16) 



(17) 



where Eh is the barrier height and E„ is the energy of the 
n'th level. Above the surface the wave function is given 
by 



Zn (z) - Z„ (z + 6z) (z > 0) 



(18) 



so that effectively the wave function has been shifted 
downwards by a small penetration depth Sz. 

Obviously Z (z) must be continuous and differentiable 
at z = 0, giving 

Szn^— An ^ Sz„Z'„ {0) = Szr^^Ai' i-yn) (19) 

This gives the probability of the neutron being inside 
the metal 









Zn (z) 



dz 



/42 



(20) 



Provided P„ is small, the calculation is self-consistent. 
We find the penetration depth and the probability 

0.01/i,P„ 2 X 10~^ (21) 



SZn 



Both increase very slowly with n, as classical intuition 
would suggest. This gives a correction to the real part of 
the energy 

SEn = PnEb - O.OlpeV (22) 

which is less than a 1% change. Potentially more impor- 
tant is the imaginary part of the potential, which gives 
rise to a lifetime 

Tn = — 1.4 X 10^ s (23) 



2 / Zn (Z) 



5V,dz 



This is almost certainly over-optimistic, since the life- 
time of higher energy UCN's is rather less than would be 
predicted by this method 0. However, it does suggest 
that the confinement time for these states is likely to be 
at least comparable to the neutron lifetime of 885 s. 



V. MAGNETIC DIPOLE INTERACTION AND 
THE EDM MEASUREMENT 

To utilize this apparatus, we assume that neutrons 
are trapped in this lowest state with all the higher en- 
ergy states centrifuged out. The first observation is 
that it would be very easy to produce polarized neutron 
states by applying a static magnetic field. The mag- 
netic and gravitational forces on the neutron are equal 



for d ^ 

iiJi "mag 



mg or a modest ^ = l.lTm ^. Hence 



by applying such a field, one spin state would have its 
binding energy reduced, its bounce height increased and 
hence would be centrifuged out , while the other spin 
state would be more strongly bound. Such a field could 
easily be achieved by a current loop surrounding the ap- 
paratus, although the spatial variation slightly compli- 
cates matters. 

To study the state in more detail, we consider a tech- 
nique of resonantly exciting a neutron in the n= 1 state 
to the n = 2 or higher. This could be done via a physical 
oscillation of the form discussed in 0, or an oscillating 
magnetic field. 



_B(z, r) = Bq (z, r) cos {cut) 



(24) 



If the field is spatially uniform, the excitation matrix 
elements will vanish identically. We therefore assume 
that the magnetic field is produced by a loop of wire as 
above. This gives rise to an interaction matrix element 
of the form 



Tn 



{n\fI.B\l) 



(25) 



The probability for resonant transition to the first ex- 
cited state (we follow ^] is given by 



P = ^ sin^ (7t) 

7 



(26) 



where 



rji2 

-f^^n^ + suj^, = -TT,Suj^ = (uji - uj-,f - (27) 



The frequency 7 must satisfy 



7 << \uji - W2I = 254iJz 



(28) 



which implies a maximum magnetic field of a few mil- 
ligauss. Since the excited states will have an extremely 
small natural width, one could, in principle, measure the 
transition energy quite accurately. 

To make an EDM measurement requires a field as large 
as possible, since the dipole moment is very small: we 
assume i?max = 10^Vm~^ . To show what limits one 
could expect to put on the EDM, we suppose that the 
oscillating field E{z, r) = Eq (z, r) cos [tot) is supphed 
by a small charged sphere on the axis of symmetry of the 
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apparatus, at a depth zq. This gives rise to a perturbing 
potential of the form 



dei.E 



rfei-Emax jz + Zp) {Rj + zg) 
7 I o\3/2 

((r + i?i)V(z + zo)'j 



(29) 



If Jcj is set equal to zero, we obtain the transition prob- 
abilities shown in Fig [S] for 4 values of the dipole mo- 
ment dei, where we have included the neutron lifetime. 
The transition time I/7 for dez = 3 x 10"^"'^ is approxi- 
mately 1300 s, which is comparable to the neutron half- 
life. Clearly an optimized design would improve the re- 
sult, but probably not by more than a factor of 10 at 
best. Hence this technique cannot give a result which is 
competitive with the best current methods. 



VI. SHORT DISTANCE BEHAVIOUR OF 
GRAVITY 

Several models have suggested that the short distance 
behaviour of gravity may be modified from Newton's law. 
It has been pointed out by Bertolami and Nunes 9] that, 
in principle, the system provides a method of measuring 



corrections to the gravitational interaction with a range 
of a few microns. We consider whether this could be done 
in practice. The gravitational potential is taken to be 



V{r) = 



GmiTO2 



(30) 



A slab of material of density p under the apparatus 
gives rise to a extra interaction the form 



-\z 



(31) 



This, of course, depends on our assumption of the slab 
being infinite: in practice, if the thickness of the slab is t, 
the approximation breaks down for A t^^. This modi- 
fies the energy spectrum and hence in principle could be 
measured. For A << h^^ (i.e. A ^ 10^^^^ or less), this 
term simplifies to give a constant term (which can be ig- 
nored, since it shifts all the levels equally) and a term of 
the form 



SV{z) = —ttK pniNZ^ 



(32) 



which is independent of A. 

In Fig , we show the shift in the transition energy 
for the 5 lowest levels as a function of A, for K = 1. We 
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conclude that if the extra couphng constant was large 
and energy shifts could be measured down to 0.1 peV, it 
would be possible to look for an extra gravitational shift 
with a range of ~ 1^. 

One can again use resonant excitation of the ground 
state to put a better limit on the extra interaction. We 
assume that the slab is vibrated vertically with a fre- 
quency u with an amplitude zq so that the distance from 
the top of the slab to the neutron is given by 



n=l state, it will then be quite difficult to excite it and 
detect it on its escape from the centrifuge. 



z (t) = z + zq{1 ~ sin {tut)) 



(33) 



This means that the neutron would feel a time-varying 
potential of the form 



(34) 



which is obviously highly anharmonic. We can use the 
expansions [To| 



oo 

gA.o cos(^t) ^ + 2 ^ /fc (Azo) COS (ujt) (35) 

A:=0 

Since the transition energy to the n'th excited bounce 
state uji — uj„ is not a simple multiple of w = \uji — 0^2 1 
, the resonance condition means that we can ignore all 
terms except one in the expansion. This gives rise to an 
interaction of the form 



SViz,t)= sin(^t)^^^?^e-^(^+-) 



(36) 



X ((A (2z + zo) + 4) h (Azo) - 2Azo/i (Azo)) 

This gives rise to an transition time which depends 
linearly on K and non-linearly on A. 

In Fig [7\we show this time as a function of A for K=l. 
It would be possible to use this to investigate short range 
changes to the gravitational interaction, but this does not 
allow a measurement to be extended down into the the 
phenomologically interesting range of K G. 



VII. CONCLUSIONS 

We conclude that it is possible to isolate this novel 
quantum state, which, of course, has its own intrinsic in- 
terest, since the number of simple quantum mechanical 
states is small, and this, along with Rydberg states, are 
unique in having a macroscopic size. It is obviously disap- 
pointing that it does not seem possible to set new limits 
on the neutron EDM, but it does appear possible to test 
the gravitational interaction at very short distances. 

We do not underestimate the problems of performing 
a real experimental measurement. It would obviously be 
very hard to observe the bouncing neutron at all, given 
the low flux rates combined with the huge flux of more 
energetic neutrons. Even if a neutron is captured in the 




FIG. 7: Plot of the transition time t as a function of A from 
the ground state to n=2,3 and 4. The neutron half-hfe is also 
shown. 



There is no reason in principle why this method can- 
not be extended to atoms. Atomic fountain experiments 
provide atoms which can have temperatures as low as 
1 nK, and hence would provided a much more numerous 
source. The major unknown is that atom-surface Van der 
Waals interactions are likely to be fairly strong and at- 
tractive. This would seriously distort the lowest bounce 
states. However, atoms are stable and cannot easily be 
absorbed. 

Finally we note that, with the addition of an appropri- 
ate static magnetic field, it would be possible to produce 
quantum states with an even lower energy in a controlled 
fashion. 
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APPENDIX A: NUMERICAL SOLUTION OF 
THE SCHRODINGER EQUATION 

The unperturbed wave function can be written 

«'fcin(r,z,0) = /,fe(r)Z„(z)e*''^ (Al) 

The radial function has to vanish at r — Ri and r = 
R3. The first gives the solution 

fir) - Yi (Ki?i) Ji (nr) - Ji (kRi) Yi {nr) (A2) 

and satisfying the second condition defines the eigen- 
values Kki via 

Yi (KkiRi) Ji {KkiRs) - Ji (KkiRi) Yi (KkiRs) = (A3) 

This is solved numerically: it is easy to get an asymp- 
totic form for the eigenvalues: 

= -f, ^ (A4) 

We must now find the perturbed wave functions 
^m{f, z) for the problem with a step in terms of the un- 
perturbed wave- functions '^k\n{r,z) (the dependence on 
(j) is irrelevant). To satisfy the boundary conditions with 
the step, we write 



{ipk\n (r, z) -t (r, z, Zs (r)) V'fcin (f, (r))) 



(A5) 



(A6) 



or 



$™(r, z)^Y. «fcin^fc In {r, z, Zs (r)) (A7) 



The function t(r, z,Zs(r)) must satisfy 
i {r, Zs (r) , Zs (r)) = 1. so that V'fe,™ (r, Zs (r) , Zs (r)) = 0, 
must vanish as z — > ck) at least as fast as the Airy 
function itself, and must vanish as r —^ Ri. It must also 
be smooth. Apart from this it is quite arbitrary: we use 



-^^^ z<h,r<R2 



t{r,z,Zs (r)) 



where 



{R2-R1) 

(2i^ple-P'{^-'^f) z>h,r<R,;, 
-P\^-'-r),>h,r>R, 



e 



(A8) 



[iz^hf + {R,~rf) 




1/2 
1/2 



Uh = {R2 - Ri 



which satisfies these properties. 
These states have energy given by 



(A9) 



Em — 



\ m\ rai 



(AlO) 



or 



^™ E E "SaS / i^ki {r)^Pki ir)df 

kl k'l' •' 



(All) 



The principal problem is caused by evaluating 

^'^'ipk,n [r, z, Zs (r)) 

. This must in general be done numerically. We regard 
the flj, as variational parameters and use them to mini- 
mize E„i- This gives 

N ^ N 

EmY^a^Gu^ ^Y."'k{Hki+Hiu) (A12) 


where 

Gki = / tpkipidrdz, Hki = I tljkHijji{x)drdz (A13) 



which is an eigenvalue equation for the Em- Note that 
Gki is symmetric but Hki is not. However, the result- 
ing eigenvalue equation contains only symmetric matrices 
and can be solved by the techniques of pd| . 



